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Abstract 

The methods used to analyze the aerodynamic 
performance of V/STOL inlets at the NASA Lewis 
Research Center are briefly described. Recent 
extensions and applications of the method are em- 
phasized. They include the specification of the 
Kutta condition for a slotted inlet, the calcula- 
tion of suction and tangential blowing for bound- 
ary layer control, and the analysis of auxiliary 
inlet geometries. A comparison is made with ex- 
periment for the slotted inlet and also for tan- 
gential blowing. Finally, an optimum inlet diffu- 
ser velocity distribution is developed. 

Noiiienclature 

Cf skin friction coefficient 

Cp pressure coefficient 

D fan diameter 

L inlet length (0.3048 m) 

M Mach number 

m mass flow rate 

P total pressure 

5 surface distance 

S.P. stagnation point 

V velocity 

X x-coordinate 

^ y-coordinate 

X nondimensional distance from the 

start of diffusion normalized 
by the total diffusion length 

0 angle o*‘ attack 

6 inlet y. v angle 

1 displacement thickness 

9 circumferential angle (9 ■= 0° at 

windward plane) 


Subscripts 


B 

blowing plenum 

c 

control station 

de 

diffuser exit 

e 

edge of the boundary layer 

i nlet 

inlet 

J 

blowing jet 

max 

maximum 

ref 

reference 

suction 

suction 

T 

throat 

00 

free stream 

PAN 

fan face 


I nlroduct ion 


In recent years, nifny different airframe/en- 
gine configurations have been proposed for V/STOL 
aircraft. Some of the proposed configurations 
impose rather severe flow conditions on the pro- 
pulsion system inlet. For example, the approach 
ana takeoff flight paths of a tilt nacelle V/STOL 
aircraft may result in inlet angles of attack up 
to IZO'. A major concern for the inlet designer 
at these conditions is the possibility of inlet 
internal flow separation. Separation free flow is 


desired to minimize thrust loss, minimize fan 
blade stress, and prevent engine stall. To assist 
the designer in providing inlet designs with sei a- 
ration-free flow, reliable theoretical methods of 
inlet flow analysis are desired both for the in- 
let design process and to interpret and augment 
the results of wind tunnel testing. The methods 
should be able to calculate the potential and 
boundary layer flows in inlets of arbitrary geome- 
try and at arbitrary flow conditions. 

Such methods of anclysis have been developed 
over the past several years at the NASA Lewis Re- 
search Center. They consist of a series of compu- 
ter programs documented in Refs. 1 to 6. Compari- 
sons with experimental results are presented in 
Refs. 7 to 12. Since these reports have been pub- 
lished, the programs have been extended and ap- 
plied to mere difficult inlet problems. 

The present paper will briefly describe the 
basic method of analysis. The major emphasis, 
however, will be on presenting the recent exten- 
sions and applications. The topics covered in 
this paper are; the flow about a slotted inlet; 
the flow in an auxiliary inlet; the analysis of 
suction and blowing boundary layer control; and 
the development of an optimum diffuser velocity 
distribution. 

Basic Method of Analysis 

The basic problem to be solved is to calculate 
the compressible viscous f 'ow in inlets of arbi- 
t> ary geometry and at arbitrary operating condi- 
tions. A series of computer programs developed at 
the NASA Lewis Research Center are used to solve 
this problem. A flow chart depicting the sequence 
for using these programs is presented in Fig. 1 
with the basic programs on the left, and recent 
extensions on the right. All programs start with 
the geometry program, upper left-hand block, which 
creates the discrete control points for each geo- 
metric configuration. Then the incompressible 
potential flow program is used to calculate the 
basic solutions to the problem. These basic solu- 
tions are combined into a solution that satisfies 
the inlet operating conditions of freestream ve- 
locity, angle of attack, and inlet mass flow. 

Next, the imcomoressible flow is corrected for 
compressibility effects. Ihe compressible poten- 
tial flow solution is then used as an input to the 
boundary layer program which calculates the lami- 
nar, transition and turbulent boundary layer char- 
acteristics, and predicts flow separation. 

Two iteration loops are available as shown to 
the left in Fig. 1. The first adds the displace- 
ment thickness to the geometry to improve the ac- 
curacy of the potential flow and boundary layer 
calculations. The second incorporates an automa- 
tic angle of attack sweep to find the separation 
angle of attack of an inlet at any given value of 
freestream velocity and inlet mass flow in one 
uninterrupted computer run. 


The recent extensions to these programs, which 
will be the major emphasis of this paper are: (a) 
to calculate the flow in an inlet with a leading 
edge slot; (b) to calculate the performance of 
suction and tangential blowing boundary layer con- 
trol concepts, and (c) to analyze the flow in aux- 
iliary inlet geometries. The method has also been 
applied to finding the optimum velocity distribu- 
tion in a subsonic diffuser. 

Geometry Program 

A program called SCIRCL is used for 2-D and 
axisymmetric geometries. For an axisymmetric in- 
let case, the geometry is represented by its meri- 
dional profile which is shown in Fig. 2(a). Both 
the external and internal ducts are extended far 
downstream so that accurate potential flow solu- 
tion can be obtained in the region of interest. 
SCIRCL breaks the profile into segments with a 
control point on each segment which are used for 
potential flow calculations. The program also 
calculates information such as curvature, wall 
angles, and flow area distribution which are very 
useful in preliminary screening of proposed inlet 
shapes. In addition to the surface points, SCIRCL 
generates off body points (like flow measuring 
rakes) also shown in Fig, 2(a) at axial locations 
where the velocity Piufile or the streamlines are 
desired. 

The 3-D geometry program, applicable to inlet 
geometries like that shown in Fig. 2(b) and dis- 
cussed in Ref. 13, allows the user to input a rel- 
atively small number of points to define the inlet 
and centerbody. The routine then enriches the 
point number and redistributes the points for good 
potential flow analysis. The detailed description 
of this geometry package is given in Ref. 5. 

Incompressible Potential Flow Basic Solutions 

The Douglas Neumann Program(^ 14-16) 
for calculating the incompressible potential flow 
field. The following basic solutions are obtaineo 
by this program: 

1. Static solution (V„ = 0) 

2. Uniform axial flow solution 

3. 90° angle of attack solution 

4. 90° angle of yaw solution (for 3-D geome- 
try only) 

In general, to obtain the basic solutions, the sur- 
face is replaced by a numiier of panels on which 
there is a surface source (or sink) distribution 
of unknown strength. For 2-D and axisymmetric 
cases, the source density can be a constant, linear 
or parabolic. For the 3-D case, only a constant 
source density can be used. The strength of source 
distribution varies over the surface in a manner 
such that at every control point the normal veloc- 
ity is zero. The best static solution is found to 
result from using a vorticity distribution on the 
cowl surfaces (as opposed to a source or sink dis- 
tribution). 

Linear Combination and Corrections 

The basic solutions obtained from the incom- 
pressible potential flow calculation are combined 
linearly into a solution of interest having arbi- 
trary flow conditions of free stream velocity, 
mass flow rate, and angle of attack. (1^) In 
cases where a Kutta conoition is required, the 
constants for linear combination are readjusted to 
satisfy the kutta condition. The linearly com- 


bined incompressible solution is then corrected 
for compressibility. ''O) if local velocity 

is supersonic, it is further corrected by an em- 
pirical supersonic correction formula. I'9) The 
final potential flow solution can now be used as 
an input to the boundary layer program. 

Boundary Layer 

TheaTia lysis of the boundary layer uses a 2- 
dimensiorial compressible boundary layer program. 

The complete documentation of the boundary layer 
program is given in Ref. 6. The program calcu- 
lates important boundary layer parameters such as 
displacement thickness, momentum thickness, and 
skin friction coefficient, Cf. It also provides 
the boundary layer velocity profiles at any de- 
sired station. The location of transition from 
laminar to turbulent flow can either be predicted 
by the program or can be specified by the user. 

Flow separation is defined to occur when the skin 
friction coefficient becomes zero. 

Recent Extensions 

The discussion thus far has described the ba- 
sic method of analysis. Now the discussion will 
focus on describing the recent extensions which 
were motivated in part by the following thoughts. 

It is desirable to design a V/STOL inlet as short 
and as thin as possible in order to reduce the 
weight, to reduce the friction drag at cruise, and 
improve pilot visibility. However, when an inlet 
is too thin the peak velocity is so high that the 
subsequent adverse pressure gradient causes the 
flow to separate at the lip resulting in a low 
pressure recovery and high distortion. Several 
ways to help control this possible separation are 
by the use of an inlet lip slot, the use of auxil- 
iary inlets, by suction or blowing boundary layer 
control, or by optimizing the surface pressure 
distribution. The analysis techniques to analyze 
these possibilities are considered next, starting 
with the slotted inlet. 

Slotted Inlet 

An axisymmetric slotted inlet is shown in 
Fig. 3. Two cases are considered, zero degrees 
angle of attack and angle of attack. At zero deg- 
rees angle of attack or at zero forward velocity 
for an axisymmetric geometry the flow is axisym- 
metric and specifying the kutta condition (i.e., 
the flows from the upper surface ana the lower 
surface join smoothly at the trailing edge.) at 
one circumferential location is equivalent to spe- 
cifying the Kutta condition around the entire cir- 
cumference of the slat. Calculations were maoe 
for static conditions, V» = 0, and the results 
are shown in Fig. 3. Experimental data arc inclu- 
ded for comparison. The agreement is quite good 
on the main inlet cowl surface (points 6 to 9). 

The agreement between the theoretical and experi- 
mental surface velocities is also good on the slat 
(points 1 to &). Figure 3 shows that the peak vel- 
ocities occurs at point 2 on the slat and point 7 
on the main cowl. Both peaks are consioerably 
lower than the peak velocity of the inlet without 
the slot,l^^) which is also indicated on tne 
figure. Thus the addition of the slot has unloa- 
ded the lip of the thin inlet. 

The case of a slotted inlet at an angle of 
attack is more difficult, because with the present 
method, the Kutta condition cannot correctly be 



Inposed at all circumferential positions simulta- 
neously. The Kutta condition can be imposesd lo- 
cally at one circumferential location, with the 
inlet at angle-of-attack, by adjusting the mass 
flow rate through the slot. Results obtained us- 
ing this approach are shown in Fig. 4. The agree- 
ment between theory and experiment at these condi- 
tions is quite good, however, it should be noted 
that the comparison is being made at the 270° 
circumferential position (side of the inlet). At 
this particular location, angle-of-attack effects 
tend to be minimal and in fact the surface veloc- 
ity distribution does not change much at all as 
angle-of-attack is increased from 0°. At other 
circumferentiai locations, a local specification 
of the Kutta condition has not been a very suc- 
cessful analysis method and work is now underway 
to extend the procedures to permit specification 
of the Kutta condition around the entire circum- 
ference of the slat at any angle-of-attack. Ex- 
tending the method of analysis to include this 
feature will result in the ability to analyze a 
new class of inlet geometries, specifically those 
that employ leading-edge slats and slots. 

Auxiliary Inlet 

The method of analysis has also been extended 
to include auxiliary inlet geometries. Auxiliary 
inlets increase the total inlet flow area thereby 
reducing the amount of airflow that must be taken 
into the main inlet. It is another technique for 
preventing flow separation on cowl lips at static 
and low flight speed conditions. 

An auxiliary inlet on the top of a conventional 
subsonic inlet and an auxiliary inlet on the top 
of a supersonic inlet are considered. A conven- 
tional subsonic inlet with an auxiliary inlet is 
shown in Fig. 5. A continuous "N-line“ (in the 
longitudinal direction) is required for the cur- 
ren*: version of the 3-D potential flow program. 

When the N-line meets the auxiliary inlet opening, 
it is rerouted along the side wa'l of the auxil- 
iary inlet and then proceeds back to the original 
N-line as shown in Figs. 5(a) and (b). Additional 
N-lines are added to completely panel the inlet. 
This particular example required 682 panels to 
describe the geometry. 

Figure 5 also illustrates the technique used 
to yield the best static solution for an inlet 
with an auxiliary inlet: Two inlet-duct systems 

are considered, one with a straight duct and one 
with a flared duct as shown in Figs. 5(a) and (b). 
The flared duct induces more flow through the main 
and auxiliary inlets. The difference between the 
velocities for the flared inlet and nonflared duct 
for a free stream uniform flow then provides the 
static solution. This procedure was adopted be- 
cause the velocities in the region of an auxiliary 
inlet were unreal istically large when the vorti- 
city distribution method noted earlier, was used 
for the static solution. 

The computer tir.ie for the basic solutions with 
the 682 panels is quite high, 19 minutes. How- 
ever, the basic solutions are only computed once 
and are stored in the computer for later use in 
obtaining solutions of Interest. Subsequent cal- 
culations using a linear combination method re- 
quired only 5 seconds of computer time. 

In order to show the effectiveness of an aux- 
iliary inlet for reducing the inlet ak velocity. 


a calculation was made for a conventional inlet 
with and without an auxiliary inlet as shown in 
Fig. 6. The peak velocity ratio is reduced con- 
siderably from 3.4 to 2.5 at the highlight of the 
windward plane when the auxiliary inlet is open 
which is, of course, the desired result. 

The pressure coefficient distribution on the 
surface around and within an auxiliary inlet is 
illustrated in Fig. 7 . In this particular case, 
the results are shown for a supersonic inlet with 
an auxiliary inlet on the top as illustrated in 
Fig. 7(a). (This i igure also illustrates the hid- 
den line removal feature of the geometry package.) 
For this case, the Kutta condition is satisfied at 
the trailing edge of the front surface (point A), 
and the high velocity occurs at the lower portion 
of the downstream surface of the auxiliary inlet 
(point B). By examining pressure distributions of 
this type, this analysis method can be used to 
study such variables as the shape, the size, num- 
ber of, and location of the auxiliary inlets. 

Suction and Tangential Blowing 

Another recent extension to the basic methods 
is the analysis of suction and tangential blowing 
boundary layer control systems. Suction controls 
the boundary layer by removing that portion of it 
not having sufficient momentum to negotiate the 
subsequent adverse pressure gradient. Blowing 
controls the boundary layer by reenergizing it 
with a thin jet of high velocity air injected tan- 
gentially into the boundary layer. 

Some results from this analytical method are 
shown in Figs. 8 and 9 for suction and blowing, re- 
spectively. An axisymmetric inlet having a diffu- 
ser exit diameter of 0.508 m was analyzed at a free 
stream Mach number, (4. = 0.12, throat Mach number, 
M-|- = 0.4, and angle of attack, o ■= 60°. The skin 
friction coefficient distribution on the internal 
surface of the windward cowl is shown along with 
boundary layer velocity profiles at several loca- 
tions. Without boundary layer control, the solid 
line, the flow separates at S/L = 0.81 '•■here the 
skin friction coefficient becomes zero. The bound- 
ary layer profile, just before separation, is quite 
weak compared to the one upstream at S/L - 0.48. 
Separation is prevented when the boundary layer is 
controlled by suction (Fig. 8) - the dashed line. 

It was necessary to bleed off only 0.12 percent of 
the inlet mass flow to prevent separation as indi- 
cated by the nonzero skin friction coefficients. 

The static-to-total pressure ratio at the suction 
location is 0.796. 

For the blowing boundary layer control (Fig. 9) 
a blowing velocity ratio, jet velocity to boundary 
layer edge velocity, Vj/V,, = 1.75 was selected. 

For this case a blowing mass flow of 0.4 percent 
of the inlet mass flow was required to maintain 
attached flow. The reenergized boundary layer is 
clearly evident in the velocity profile just down- 
stream of the blowing slot. 

A comparison of analytical and experimental ve- 
locity profiles in the windward plane with blowing 
being used for boundary layer control was made for 
0.S08 m axisymmetric inlet at a = 70°, V]/Vq ■= 2 
and Fy/Pg = 1.2. The axial locations of the total 
pressure rakes used to measure the profiles are 
shown in Fig. 10. As indicated, the analytical 
and experimental velocity profiles agree well at 
the upstream measuring station except in the vi- 
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cinity of Y/6 - 0.375, where the analytical veloc- 
ity is too high. This discrepancy between the 
analytical and experimental results could be the 
result of the closeness of this measuring station 
to the blowing slot and the difference between the 
assumed velocity profile and the experimental ve- 
locity profile at the blowing slot exit. A similar 
trend in blowing velocity profiles was observed in 
Ref. 20. At the downstream rake location, however, 
the mixing between blowing jet and boundary layer 
has already taken nlace so that tne well developed 
analytical velocity profile agrees very well with 
the experimental velocity profile. 


optimum velocity distribution is slightly more 
consevativ',' at the beginning of the diffusion pro- 
cess than that of Stratford. Since Stratford's 
distribution is derived on the basis of zero skin 
friction throughout the diffuser, it can be consid- 
ered as a limiting case. A design velocity distri- 
bution (besides having a safety margin against 
separation) should have a slightly more gradual 
start to the pressure rise (the deceleration of 
velocity) than that of Stratford. The present 
optimum velocity distribution has this character- 
istic and deserves further investigation as a po- 
tential attractive inlet design. 


Next, a comparison of analytical and experi- 
mental separation bounds is shown in Fig, 11 for 
the same inlet. The separation parameter Vt/V« 
is used to collapse the experimental data. (Refs. 
21 and 22). Agreement between the analytical and 
the experimental separation bounds is iound to be 
excellent up to a = 80°. For angles of attack 
greater than 80°, the analytical separation p:in. 
is predicted to be ahead of blowing slot and the 
analysis can then not proceed since the boundary 
layer code cannot perform calculations beyond tne 
separation point. 

Optimum Diffuser Velocity Distribution 

Another application of the method is concerned 
with finding the optinium velocity distribution in 
a suDSonic diffuser. This velocity distribution 
wi''! result in the shortest no-bounda'*y layer con- 
trol inlet and the lowest loss for the required 
amount of diffusion. 


Concudinq Remarks 

An analysis method based on incompressible 
potential flow corrected for compressibility was 
described. Several sample calculations compared 
well with experimental data. The most recent ap- 
piic'tions include an inlet with a leading edge 
slat, an auxiliary inlet, and suction or blowing 
boundary layer control. Ar optimum diffuser velo- 
city distribution was also developed. This paper 
shows that the present methods can be a very power- 
ful tool for the analysis and the design of V/STOL 
inlets. 
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GEOMETRY PACKAGE 
(2D, AXI, 3D) 



Figure 1. - Schematic diagram of computer programs. 
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